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The structural weight reduction of a pressurized Hybrid Wing-Body (HWB) fuselage is a serious challenge. 
Hence, research and development are presently being continued at NASA under the Environmentally 
Responsible Aviation (ERA) and Subsonic Fixed Wing (SFW) projects in collaboration with the Boeing 
Company, Huntington Beach and Air Force Research Laboratory (AFRL). In this paper, a structural 
analysis of the HWB fuselage and bulkhead panels is presented, with the objectives of design improvement 
and structural weight reduction. First, orthotropic plate theories for sizing, and equivalent plate analysis with 
appropriate simplification are considered. Then parametric finite -element analysis of a fuselage section and 
bulkhead are conducted using advanced stitched composite structural concepts, which are presently being 
developed at Boeing for pressurized HWB flight vehicles. With this advanced stiffened-shell design, 
structural weights are computed and compared to the thick sandwich, vaulted -ribbed-shell, and multi- 
bubble stiffened-shell structural concepts that had been studied previously. The analytical and numerical 
results are discussed to assess the overall weight/strength advantages. 
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Nomenclature 

spacing between y-directional stringer or frame. 

beam column cross section area. 

section area of single v-directional frame or stringer. 

area of single y-directional stringer or frame. 

spacing between v-directional frame or stringer. 

total breadth of stiffened panel between end supports. 

bending rigidity of plate £^ 3 /12(l - i?). 

bending rigidity of stiffened orthotropic plate about the y axis. 

bending rigidity of stiffened orthotropic plate about the x axis. 

Young’s modulus of orthotropic plate in x and y directions, 
yield stress in compression along x and y directions, 
yield stress in tension along x and y directions, 
shear modulus. 

VD x D y . 

area moment of inertia for beam bending about neutral axis, 
area moment of inertia of x- and y-stiffeners about neutral axis, 
total length of stiffened panel between end supports. 

Possion’s ratio along x and y axis. 

running bending moments about x and y axis, respectively (lb-inch/inch). 

running in-plane load along x and y directions, respectively (lb. /inch). 

compression load on beam-column. 

compression buckling load on beam-column. 

cabin pressure of 9.2 psi (2 P =18.4 psi). 

running normal load on beam or pressure on plate. 

shell or skin thickness. 

beam or plate deflection at mid-point. 

v and y reference axes of the plate, respectively. 

normal distance from beam or plate neutral plane. 

neutral axis location of x- and y-stiffeners from skin mid plane. 
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I. Introduction 


D 13 

J -^esign of efficient pressurized fuselage configurations for the Hybrid- Wing-Body (HWB) and Advanced 

Mobility Concept 4 (AMC-X) flight vehicles are challenging opportunities for the efficient structures technology 
development. The stiffened cylindrical cross-section fuselage of Boeing 111 and 787, shown in figure 1, are quite 
efficient, since, when pressurized, the skin is mostly under tensile membrane stress. The Airbus A3 80 has an elliptic 
egg shaped cross-section which is strengthened by the twin deck floor beams. However, in a HWB fuselage with 
rectangular cross-section, the internal cabin pressure combined with flight loads primarily result in highly non-linear 
bending stress and large deformation. It can be shown 5 ' 6 that these bending stresses are theoretically one order of 
magnitude higher than the membrane -hoop stress in a conventional cylindrical fuselage of similar size and skin 
thickness. Moreover, resulting deformation of the outer aerodynamic surface could significantly affect performance 
advantages provided by the large lifting body. Thus, it is necessary to design an efficient fuselage structure in order 
to reduce the overall structural deflection and weight penalty, while satisfying the design stress, strain, and buckling 
safety factors, under the critical design loads 7 " 10 . Many structural concepts such as the conventional skin -stringer- 
frame, ribbed-double-shell, multi-bubble stiffened shell, thick sandwich, J-Frame, Beaded-Hat, and stitched 
sandwich frame composite construction have been studied under AFRL and NASA research projects 8 " 10 . These 
concepts are schematically shown in Figure 1 . 
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Figure 1. Fuselage section concept analysis models of B777, A380 and HWB class of vehicles. 


A highly innovative advanced technology acronymed PRSEUS (Pultruded Rod Stitched Efficient Unitized 
Structure) is being developed at Boeing Phantom Works, with application to the pressurized HWB vehicles 11 ' 16 . 
This involved analysis, design, fabrication and testing of PRSEUS coupons, panels and fuselage components, under 
the NASA Environmentally Responsible Aviation (ERA) Program’s Airframe Technology subtask for structural 
weight and risk reduction. This development process 1718 for the PRSEUS concept is schematically shown in Figure 
2. Some initial results at the coupon, panel and fuselage section under one critical overpressure load condition were 
presented in Ref. 18. This paper will discuss additional analytical and numerical results of parametric design studies 
of conceptual fuselage section and bulkhead for comparative weight analysis and optimization. Alternate concepts 
are also considered for systems study in an attempt to improve the present design farther. Only local idealized 
fuselage section is considered in this paper, with a limited set of critical design loads. The objective is to minimize 
the local structural weight by 10-20% with stitched composite construction combined with novel structural 
configuration, compared to the sandwich and multi-bubble concepts that were studied previously. At the vehicle 
level, aerodynamic and structural systems study for the 100, 200, 300, and 400-passenger versions were considered 
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in Refs. 19-23, in order to estimate the full vehicle structural weight at all critical load conditions. For a recent multi- 
bubble HWB configuration study, see Ref. 24. 
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Figure 2. PRSEUS concept development cycle for analysis, design and testing. 


II. Classical Stiffened Plate Analysis 

Since the Finite Element (FE) model development of the PRSEUS concept requires many levels of simplification 
and idealization for extension from the coupon level to the fuselage level 18 , an analytical model is developed at the 
panel level for rapid parametric study and sizing. This would also enable one to evaluate the effects of skin 
thickness, frame and rod-stringer spacing and section geometry on the panel deflection due to a pressure load and for 
buckling analysis due to in-plane compression load. Classical plate theory 25 ' 26 and equivalent orthotropic stiffened 
plate equations 27 ' 29 were applied to simplify the PRSEUS model analysis for a rectangular plate. These simple 
spread-sheet based analytical results are also used for sanity check of the high fidelity FE analysis of panels, 
fuselage section and bulkhead. SolidWorks 30 was used for rapidly building detailed FE models of the vehicle 
fuselage sections for structural analysis and design. A large number of design options were considered such as the 
orthotropic skin thickness, frame and rod-stringer spacing, for obtaining stress and strain distributions, deformation 
pattern and for identifying critical areas. For each configuration, structural weight was computed and compared. 



Stiffened Plate Bending Analysis: Consider the biaxially stiffened plate shown in Fig. 3. For an isotropic thin 
rectangular plate of uniform thickness t, length L and width B, the maximum deflection can be expressed as in Eq. 1 
which is a function of the pressure load q, plate bending stiffness D and the factor k that is calculated from analytical 


3 





solution for this boundary condition. The factor k is presented in tabular form 25 , as a function of the ratio B/L 
(B>L). For example, for a square plate with built-in edges, k=0. 00126 for B/L=l (Table 1). For a simply supported 
isotropic square plate, k = 0.00406. 

(1) w max =k.qL 4 /D where D = Et 3 /12(l-\?) 

Let us assume that the uniform plate-theory assumptions are also applicable to this equivalent stiffened plate. From 
ideal orthotropic plate theory 25 , for a simply supported stiffened rectangular plate, the closed form solution for 
deflection w given by Eq. (2). 
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The running stifnesses D x , D y and H can be approximately defined by Eqs. (3-5). 
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(5) H = jD x D y 


The terms in Eq.(3-4) are explained in Figure 3. The equivalent plate bending stiffnesses D x , D y resists bending 
moment M x , M y about the v-and y- axes, respectively. The attached longitudinal stiffening frames, each with area A sx 
and individual bending stiffness E X I X move up the cross-sectional-area center by Z ox from the plate neutral (mid) 
plane, thereby increasing the total area and the bending stiffness. The symbol E x is used to differentiate the stiffener 
modulus from the plate modulus E x . Similarly, the transverse stringers or rod-stiffeners, each with area A sy and 
bending stiffness E y I y , increase the bending stiffness D y by moving up the transverse-section-area center by Z oy from 
the plate neutral plane. 

For the stitched composite orthotropic skin, the following material properties were used 8-9 , namely E x = 9.75x1 0 6 
psi, E y = 4.86xl0 6 psi, v x = 0.39, and v y = 0.2. A stitched composite frame with a O.linch skin wrapped around a 6 
inch high, 0.5 inch thick foam core and a 4-in. wide flange and a 0.1 inch thick base cover strap provide a cross- 
sectional area of approximately 2 in. 2 , without considering the foam core, that is mostly for shape stabilization. The 
frame dimensions are schematically shown in Figure 4(a). The frame flange is stitched to the base skin and are 
spaced at b=20 in. interval, and with a 0.1 in. 
skin thickness, the skin area between the 
frames is approximately 2 in 2 . Together, the 
frame and the skin segment between 
adjacent frames have running bending 
stiffness D x of approximately 5.8xl0 6 lb-in. 

By using a 0.2 inch thick frame cap (Fig. 

4b), D x increases to 6.3x1 0 6 lb-in. This 
would reduce the stresses, which are highest 
in the frame cap. For the first rod-stringer 
shown in Fig.4c, with base to rod-center 
height of 1.4 in., D y = 0.75x1 0 6 lb-in. 
including a 2 in. cover-strap and base skin 
between the rod-stringers spaced at 6 inch 
apart. For the rod-stringer in Fig.4d, with 
base to rod-center height of 1.65 in., 

D y =1.08xl0 6 lb-in. The rod-stringer in 
Fig.4e, with 0.208 in. flange thickness, and 
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Figure 4. Frame and rod-stringer configurations. 






1 \ 
$ 


8 


1— 




1 5 


s 

§ 


8 

ui 



! V 



1 

1 8 

- 

i 



4 



base to rod-center height of 1.65 in., D y =1.45xl0 6 lb-in. Next, these frame- stiffened plate stiffness estimates were 
applied to the classical isotropic plate bending equations. The D x , D y and H that were computed for this frame- 
stiffened panel are now applied to a classical plate bending equation subject to a normal, uniform distributed 
pressure load of q psi with simple-supported and clamped edge conditions 25 ' 26 . 

Simply Supported Plate: Using Eq. (2) to an anisotropic rectangular plate with q 0 =18.4 psi , L=120 in., B=90 
in., D x = 5.77xl0 6 lb-in and D y = 0.75xl0 6 lb-in., the maximum deflection is 4.5 in. Here, stiffer frames (D x ) 
running along the longer side L, as shown in Fig. 3. If the frames are running along the shorter side B, and the rod- 
stringers are running along the longer side L, then the deflection can be obtained by simply interchange D x with D y 
in Eq.(2). The deflection at mid-point is now reduced to 2.1 inches. 

If we apply Eq. (2) to a simply supported isotropic square plate (. L=B , D x =D y =D) the maximum deflection is can 
be simply written as in Eq. (1), w mSLX =kqL 4 /D with k = 0.00406. If one takes just the first term ( m=l , n=l ) in Eq. (2), 
i.e. w max = (16/n)7 (D/L 4 +2H/L 2 B 2 +D/B 4 ), then for a square plate k-4/n -0.00416. With q=18.4 psi, D x = 5.8xl0 6 
lb-in and D y = 0.75xl0 6 lb-in, application of the exact Eq. (2) gives w max =3.23 inches, for a 104-in. square panel. 

For an orthotropic plate, one may modify Eq. (1) to w max = kqL 4 /^(DxDy) by replacing D with sl(D x D y ), with 
applicable k value. For example, with this approximation, for the same 104-in. square panel, the maximum 
deflection w max =3.52 in., about 9% higher than that from Eq. (2). This simple approximation, w=kqL 4 ND x D y may 
provide reasonable results for maximum deflection, for appropriate k value from Table l 25, 26 . 


Table 1: Clamped Isotropic Plate max deflection with uniform pressure q: 20 6 0.104 1.65 0.104 0.479 


Uniformly loaded Clamped isotropic stiffened plate: Table 35, p 202, nu 0.3 figure 91 p 197 Timoshenko 


CASE 


length 

width 

stiffness 

load 


max deflection 

shell 

NO 

max deflection equation 

L (L< B) 

B 

Dx, Dy or sq 

q (psi) 

k 

wmax=kqL A 4/D 

skin t 

1 

wmax=kql_ A 4/Dx, L=B 

120 

120 

5,805,106 

18.4 

0.00126 

0.83 

0.104 

2 

wmax=kqB A 4/Dy, L=B 

90 

90 

1,087,872 

18.4 

0.00126 

1.40 

0.104 

3 

wmax=k.q.(L A 4)/sqrt(Dx*Dy), L=B 

104 

104 

2,513,008 

18.4 

0.00126 

1.08 

0.104 

4 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B/L=1.3 

90 

117 

2,513,008 

18.4 

0.00191 

0.92 

0.104 

5 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B/L=1.33 

90 

120 

r 2,513,008 

18.4 

0.00196 

0.94 

0.104 

6 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B/L=1.5 

80 

120 

2,513,008 

18.4 

0.00220 

0.66 

0.104 

7 

wmax=kql_ A 4/Dx, B=2L B/L=2 

60 

120 

5,805,106 

18.4 

0.00254 

0.10 

0.104 

8 

wmax=kql_ A 4/Dy, B=2L B/L=2 

90 

180 

1,087,872 

18.4 

0.00254 

2.82 

0.104 

9 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B=2L 

90 

180 

2,513,008 

18.4 

0.00254 

1.22 

0.104 

10 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B = 1.55 L 

90 

140 

2,513,008 

18.4 

0.00225 

1.08 

0.104 

11 

max=4*k.q/[ Dx/L A 4 +2H/ ( L A 2. B A 2) + Dy/B A 4] 

80 

120 

5,805,106 

18.4 

0.00220 

0.80 Dx -> L (L<B) 0.104 


B/L=1.5 (L<B) switch Dx and Dy 

80 

120 

Dx and Dy 





12 

nax=4*k.q/[ Dy/L A 4 +2H/ ( L A 2. B A 2) + Dx/B A 4] 

80 

120 

1,087,872 

18.4 

0.00220 

1.48 Dy -> L (L<B) 0.104 


Table 1: Clamped isotropic plate maximum deflection with single wrap rod-stringer at 6 inch spacing and double wrap 
frames at 20 inch spacing with uniform normal pressure of 18.4 psi. 


CASE 

fsp 24 rsp 8 
t = 0.104 

fsp 24 rsp 8 
t=0.052 

fsp 24 rsp 6 
t = 0.104 

fsp 24 rsp 6 
t=0.052 

fsp 20 rsp 6 
t = 0.104 

fsp 20 rsp 6 
t=0.052 

fsp 20 rsp8 
t = 0.104 

fsp 20 rsp 8 
t=0.052 

max deflection equation 


NO 

1 

0.97 

1.07 

0.97 

1.07 

0.83 

0.92 

0.83 

0.92 

wmax=kqL A 4/Dx, L=B 


2 

1.77 

1.89 

1.40 

1.47 

1.40 

1.47 

1.77 

1.89 

wmax=kqB A 4/Dy, L=B 


3 

1.31 

1.42 

1.16 

1.26 

1.08 

1.16 

1.21 

1.32 

wmax=k.q.(L A 4)/sqrt(Dx*Dy), L=B 

4 

1.11 

1.21 

0.99 

1.07 

0.92 

0.99 

1.03 

1.12 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B/L=1.3 

5 

1.14 

1.24 

1.02 

1.10 

0.94 

1.02 

1.06 

1.15 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B/L=1.33 

6 

0.80 

0.87 

0.71 

0.77 

0.66 

0.71 

0.74 

0.81 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B/L=1.5 

7 

0.12 

0.14 

0.12 

0.14 

0.10 

0.12 

0.10 

0.12 

wmax=kql_ A 4/Dx, B=2L B/L=2 

8 

3.56 

3.81 

2.82 

2.96 

2.82 

2.96 

3.56 

3.81 

wmax=kqL A 4/Dy, B=2L B/L=2 

9 

1.48 

1.61 

1.32 

1.42 

1.22 

1.32 

1.37 

1.50 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B= 

=2L 

10 

1.31 

1.43 

1.17 

1.26 

1.08 

1.17 

1.21 

1.33 

wmax=k.q.(L A 4)/sqrt(Dx.Dy), B 

= 1.55 L 

11 

0.95 

1.05 

0.91 

1.01 

0.80 

0.89 

0.83 

0.92 

wmax=4*k.q/[Dx/L A 4 +2H/ ( L A 2. B A 2) + Dy/B A 4] 


B/L=1.5 (L<B) switch Dx and Dy 

B/L=1.5 (L<B) switch Dx and Dy 

12 

1.80 

1.96 

1.60 

1.73 

1.48 

1.61 

1.66 

1.81 

wmax=4*k.q/[Dy/L A 4 +2H/ ( L A 2. B A 2) + Dx/B A 4] 


Table 2. Parametric maximum deflections for the 12 cases in Table 1 for a combinatorial set of frame spacing 
(fsp=20 or 24 in.), rod-stiffener spacing (rsp=6 or 8 in.) and skin thickness (t=0.052 or 0.104 in.). 

Clamped Edge Plate: One can apply this approximation to the Eq. (1) for a clamped edge plates (Ref. 25, Figure 
91, pp.197, and Table 35, pp. 202). Table 1 shows the maximum deflection of typical rectangular clamped panels with 
D x =5.8xl0 6 lb-in., D y =1.08xl0 6 lb-in. (Frame in Fig.4a, rod-stringer in Fig. 4d, respectively), and 18.4 psi pressure 
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load. The equivalent plate bending stiffness D x or D y were computed from Eqs. (3-4) from classical plate theory with 
20 inch frame spacing, 6 inch rod-stringer spacing, 0.104 inch skin and frame- wrap thickness. 

The column at left shows the equation used for each of the 12 cases. Case 1 is for a frame stiffened (D x ) square 
plate in both directions x and y. Case 2 is for rod-stringer stiffened (D y ) square plate in both directions x and y. In the 
cases 3 to 6, 9 and 10, stiffness D is replaced by V(D x .D y ) as an equivalent orthotropic approximation, as discussed 
before. In case 3, for a 104 in. square plate, the maximum deflection is approximately 1.08 in. with q-18.4 psi. In case 
10, a 90x140 stiffened plate (lower bay side-cargo panels) this approximation yields a maximum deflection of 1.08 
inches. These approximations yield a median deflection value that cannot distinguish between running direction of 
frames and rod-stringers in rectangular panels. Hence, in case 11, D/L 4 in Eq. (1) is replaced by 
0.25 [D/L 4 +2H/(L 2 B 2 )+D/B 4 ] with k factors for the specific B/L ratio (L<B) from Ref. 26, for a better orthotropic 
approximation. This equation represents stiffer frames ( D x ) running along shorter direction L and the maximum 
deflection is 0.8 inches for an 80x120 panel. In case 12, D/L 4 is replaced by 0.25 [D/L 4 +2H/(L 2 B 2 )+D X /B 4 ], 
(Interchange D x with D y ) which represents stiffer frames (D x ) running along the longer direction B. In this case the 
maximum deflection is 1.48 inches. This characteristic is similar to that observed from the exact Eq. (2) for the simply 
supported orthotropic rectangular plate case. Of course, for square orthotropic plate, both equations yield same 
deflection. 

Table 2 (continuation of Table 1) shows maximum deflections for each of the 12 cases, which are computed for 
eight combinatorial design scenario of frame spacing of 20 and 24 inches, rod-stringer spacing of 6 and 8 inches and 
skin thicknesses of 0.052 and 0.104 inches. These spread sheet based calculations can be used for panel sizing, weight 
and configuration optimization as well as for sanity check of numerical finite element results. 

Ideal Beam-Column Buckling analysis: Consider an ideal straight beam-column to represent a flat cabin roof 
panel segment, with axial in-plane compression load P , and distributed running load q due to the normal cabin 
pressure. Let El be the beam bending stiffness over the span-length L. From the beam-column theory, the maximum 
deflection, bending moment, and bending stresses at mid-span are given by Eqs. (6-8). 


(6) 

Wmax = (qL 4 /384EI)T]( M ), 

where q{p)-12{2.sec{p) -2- p)/(5p) 

(7) 

Mmax (7/7/ /8) + Pc'^max 

, where // = j 

(8) 

cr = M, na x -// + P/A 



The critical bucking load P cr for a simply supported straight beam without initial deformation or eccentric loading is 

given by P cr = Then ju = ( 7i/2)V(P c /P cr ), and thus these equations are valid for P c < P cn or ju < n/2. Based 

on the beam-column equations, the stiffener frame dimensions can be sized to provide adequate section area and 
bending stifnesses to prevent buckling. In these beam bending calculations, in order to represent a stiffened panel, 
the total bending stiffness El in beam-column equation is replaced with D x b, and the total spanwise running load q is 
multiplied by b , where b is the frame spacing. Applying the simply supported beam-column Eqs. (6-8) with L=120 
in., b=20 in. D x =5.8xl0 6 lb-in. (Frame in Fig.4a), and a 9.2 Ib/in. normal running load, the maximum linear 
deflection is 4.28 in. and the maximum bending stress on the frame top is 6,000 psi without compression load. The 
critical beam-buckling compression load is 3960 lb/in., assuming that the frame and the plate segment bend like a 
beam. With a compression load P c =1000 lb/in., P/P^O.25 and the maximum nonlinear deflection is 5.72 in. The 
maximum bending plus axial stress is 49,000 psi including non-linear bending effect due to the 1000-lb compression 
load. Hence, for a simply supported ideal frame- stiffened beam-column, the buckling strength needs to be increased 
significantly. 

In the real situation, the PRSEUS panel will be cross- stiffened and all edges will be close to the clamped 
condition. The beam deflection for the clamped case can be computed by applying equal end-moments M 0 , which 
would reduce the end- slopes to zero. This moment and the corresponding maximum deflections are given by 26 Eqs. 
(9-10). 

(9) M 0 = ~(q.L 2 /12). p. x( p)/tan( p), where x{p)=3(tan{p)-p)/p 

(10) w max = ( MoL 2 /8.E1).2( 1 -cos( p) )/ju cos( p) 

For the clamped end condition, the linear buckling load P cr =4.7rEI/L 2 or 15,000 Ib./in. in-plane compression running 
load. The end moment M 0 is calculated to be 230,000 lb-in , and the corresponding maximum bending stress at the 
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frame-top is approximately 78,000 psi at beam ends. The maximum non-linear deflection at mid-span is then 4.8 
inches. In order to reduce the deflection and alleviate this high local stress, the panel skin and frame ends need to be 
reinforced significantly. 


Panel Buckling Analysis: Consider the simply supported stiffened plate in Fig. 3 under normal loading q and in- 
plane compressive running load N x , then the solution for deformation can be expressed as a modified form of Eq.(2) 

, N 16 q n ^ ^ sin( m7Dc / L)sin(n7ry / B) 

(ID w(x,y) = — f- 2, 2.- 


7T 


m=l,3,5 n=l,3,5 m 4 „ 2 mV „ ft 4 „ N 


mn (-^ D * + 


2m V H + nA D 


2 r 2 

71 L 


) 


The value of N x that renders the denominator of Eq. (11) to zero will provide the in-plane buckling load. Thus for 
m=l, n= 1 , in-plane compressive running load N xcr = n L [D/L + 2*H/(L .B ) + D/B ]. For the orthotropic 
equivalent plate with L=120 in., B=90 in., D x = 5.80xl0 6 lb-in and D y = 1.08xl0 6 lb-in, the critical in-plane 
running load N xcr is 12,000 Ib/in. For an isotropic square plate, the critical minimum buckling is given by N xcr = 
4k D/L 2 . For a 104x104 in. square simply supported plate with D replaced by V(D x D y ), N xcr = 9100 Ib/in. 


III. Fuselage Section Analysis. 


Design conditions: The following critical design load conditions were studied for the fuselage and bulkhead. 

1. 18.4-psi internal static over-pressure at ground level (i.e., 2P case). 

2. 2.5-g climb condition with an internal cabin pressure of 9.3 psi (i.e., lP-2.5g case). 

3. Panel in-plane buckling load of 5000 Ib/in. with 9.3 psi normal cabin pressure. 

The maximum aerodynamic load for the lP-2.5g case also produces a maximum compression load on the 
fuselage top panels and a maximum tensile load on the fuselage bottom panels. This bending load is idealized as 
approximately 5000-lb/in. running compression and tension load along the top and bottom panels, respectively. This 
in-plane ultimate design load was derived from full-vehicle FEM analysis results 14 ' 16 and includes all safety margins. 

Material Properties: For the orthotropic stitched PRSEUS panel, the elastic modulus in the v and y directions 
are E x = 9.75xl0 6 psi, E y = 4.86xl0 6 psi, u x = 0.39, u y = 0.2, respectively. The average mass density is 0.057 lb/in 3 . 
For the foam core, E = 21000 psi, G = 7950 psi , and the mass density is 0.00361 Ib/in . The web and the flange of 
the pultruded rod-stringer have the same property as the stitched composite material. The pultruded rod has an 
elastic modulus of E'= 1.9xl0 7 psi and u = 0.29. This property is used when the built-up rod-stringer is modeled as 
a beam. In some of the previous panel analysis 18 the web was modeled as shell and the top bulb was modeled as a 
beam. The stronger x-direction of the orthotropic skin and frame-wrap are always along the frame running direction. 

Failure Criteria: In the linear analysis, the maximum allowable principal strain is 0.007 for the orthotropic skin 

and frame wraps. For 
the orthotropic skin and 
built-up pultruded rod- 
stringer average yield 
stress of F^ = 46500 
psi is used for a 
conservative safety 
factor calculation. For 
solid models of the 
foam-core, yield stress 
F y = 440 psi. The 
allowable failure 

stresses of the 
orthotropic shell are F tx 
= 105100 psi, Ffy = 
46500 psi, in tension, 
and F cx = 79200 psi and 
E cy = 37900 psi in 
compression, 
respectively. The 

maximum allowable 



Figure 5. Multibay model with 5 frames at 20 inch spacing, and rod-stringers at 8 in. 
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shear stress is 29900 psi. When the rod stringers were modeled as beams, the maximum combined bending and axial 
stresses were used to check for beam failure. The directional stresses and principal strain distributions along with 
yield stress based factor of safety were plotted to identify local high stress and strain failure areas. 

Multibay Fuselage Section: For a detailed description of the Multibay development, fabrication and combined 
loads test plan and progress, see Refs. 11-16. Figure 5 shows a multi-bay concept models along with the frame and 
rod-stringer dimensions. This model is substantially simplified, from those described in refs. 11-16, in order to 
facilitate conceptual study and what-if type analysis. This 3 bay section has a spanwise width of 360 inches, with 90 
inch high passenger bay, 72 in. high central cargo bay. This fuselage section has five foam-core frames at 20 in. 
spacing. Each frame is built-up with 6 in. high, 0.5 in thick foam-core, that are wrapped with 0.104 in. thick (2 
stack) composite skin and flanges. The flanges are stitched to the base skin and frame cover straps. The rod-stringers 
are spaced 8 in. apart on upper and lower surface of each bay, side walls, lower left and right cargo bay and 
passenger floors. Each rod-stringer is built-up with stitched composite flange, web wrapped around a pultruded 
carbon rod shown in Fig 4d. The flange is stitched to the base skin and a stringer cover strap 16 . In this multi-bay 
concept analysis, the built-up rod-stringers are modeled by beam elements with same cross-section as the built-up 
rod-stringer and cover strap. The base skin and frame wraps are modeled by shell elements. The foam core is 
modeled by solid elements. The two mid-cabin sandwich walls are modeled with 2 in. thick foam core and 0.052 in. 
stitched composite skin without cutouts. 




Figure 6. Multibay displacements (magnified by xlO) and principal strain (Five frames at 20 inch spacing, improved 
rod-stringer at 8 in. spacing, without bulkhead) under internal cabin pressure of 18.4 psi and passenger-floor 
pressure of 1 psi. 


Principal Strain El 
m 0.0119 

0.0097 


Five frames at 20 in. spacing 
with 18.4 psi pressure load, 
fixed edge BC (BAY7E) 


Displacement (inch) x 10 
1.2 


0.0076 
0.0054 
0.0033 
0.0012 
* - 0.0010 


Figure 6 shows the corresponding displacements (magnified by xlO) and principal strain under 18.4 internal cabin 
pressure on all outer walls and 1 psi pressure on passenger floors. Since the bulkhead was not modeled, an all outer- 
edge-fixed boundary condition was applied. The maximum deflection on the right and left side slanted cargo 
compartment walls is 1.2 in. Note that the cargo bays are approximated as flat panels, although in actual vehicle 
these keel panels conform to the aerodynamic shape. The directional stresses and strains on the skin are generally 
within the allowable limits, except at edge-restraints, cabin-wall-skin junctions, corner junctions and foam-core 
frame areas. The maximum principal nodal strains are of the order 0.012, at inter-cabin-junctions and unreinforced 
cabin-corners. In order to reduce large local stresses, these junctions need to be modeled more realistically with 
doubler and comer braces. These corner junction areas require diagonal braces and additional reinforcements. 
Additional design refinements are required in these critical areas. The maximum bending and axial stresses on the 
rod-stiffeners, which were modeled as beam elements are generally within allowable limits, but may exceed locally 
at the lower cargo bay due to large bending and edge constraints. The stress safety factors are locally 0.5 in this area 
on the rod-stiffeners but are well above 1 elsewhere. These high stresses can be reduced by adding corner braces, 
rod-stiffeners with a higher bending stiffness and frames with thicker frame-caps in these areas along with at least 
0.1 in. skin due to impact damage concerns. The maximum deflection values from the FE analysis are close to the 
empirical average values, shown in Tables 1 and 2 (cases 6 and 10) and generally consistent with the analytical 
deflection values of similar size stiffened plates. However, there are still some unresolved issues about the analytical 
solution of orthotropic clamped plates, including buckling analysis, which will be addressed in future. An isolated 
test bulkhead analysis is presented next. 
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IV. Bulkhead Analysis 




0.566 


0.377 


0.189 


one stack rod-stringer and Frame(not to scale) 0.0 


Two stack rod-stringer and Frame(not to scale) 


Figure 7. Bulkhead displacements under an 18.4 psi normal pressure load with fixed end boundary conditions: (a) 
one stack rod-stringer (Fig. 4d) and (b) two stack rod-stringer (Fig. 4e) with 1.65 inch base to center height. 


Fig. 7(a) shows the displacement of an isolated bulkhead configuration with 20 horizontal rods at 8 inches 
spacing, 19 vertical frames at 20 inches spacing, two stack skin and two stack frame-wrap. Under the 18.4 psi 
pressure, maximum deflections are of the order 1.13 inches on the upper bulkhead left and right bay segments. The 
principal strain is generally under 0.007 inch/inch. Local combined maximum axial and bending stresses on the 
horizontal rod-stringers that were modeled as beam are under 76000 psi. Except for the local stress at the bulkhead- 
wall joints, due to end-fixity conditions, the majority of resultant stresses are under 70,000 psi. These high stress and 
deflection can be reduced significantly, by reinforcing the rod-stringers with two-stack wrap as discussed earlier. 
This is shown in Figure 7(b). The maximum deflection is now reduced to about one inch. The outer bay width for 
this configuration is reduced by 6 in. on each side. The local maximum beam bending stresses are now under 80,000 
psi. The local maximum principal strains are now under 0.007 in/in. The deflections are close to the analytical 
values presented in Tables 1 and 2 (cases 5 and 11). However the total weight increases slightly to 878 lb. (2.59 
lb/sq. feet of surface area), as shown in Table 3. 


Displacement (inch) with 
fixed boundary condition all 
outer edges, upper & lower 
walls and cargo floor joints 


Displacement (inch) 

1.132 


Displacement (inch) with two-stack 
rod-stringer and fixed boundary 
condition all outer edges, upper and 
lower walls and cargo floor joints 


Displacement (inch) 


0.755 


0.669 


new two stack 1.65b2c rod_dim 


weight calculation Bulkhead2T (2 stack skin, 2 stack wrap, strap, 8 inch rod spacing, 24 in frame spacing) Study 2 


BULKHEAD2b 

width 

length 

area 

thickness 

vol. 

density 

weight/each weights 

group 

study 2 (SW13) 

inch 

in 

in A 2 

in 

in A 3 

1 b/i n A 3 

lbs 

X no items 

weight 

skin group 










292 

top panels 

90 

348 

31320 

0.104 

3257.3 

0.057 

185.7 

1 

185.7 


lower mid bay panel 

72 

120 

8640 

0.104 

898.6 

0.057 

51.2 

1 

51.2 


triangular panels 2 

72 

114.0 

4680 

0.104 

486.7 

0.057 

27.7 

2 

55.5 













2 stack 

rod area 

web area 

flange are 

area 

length 

density 


no of rods 


Hor rod-stringer group 

in A 2 

in A 2 

in A 2 

in A 2 

in A 3 

1 b/i n A 3 

lbs 

x no 


305 

rod_stringer 1.65 in b2C 
2.1 in base @ 8 in 

0.267865 

0.260832 

0.3276 

0.856297 

360 

0.057 

17.6 

17.33333 

304.6 













Ver Frame group 

width 

length 

area 

thickness vol 

HOwf cSW 

x no 


281 

frame core at 24 inch 

6 

162 

972 

0.5 

486.0 

0.0036 

1.7 

13.1 

23.0 


frame wrap+flanges 

16.5 

162 

2673 

0.104 

278.0 

0.057 

15.8 

13.1 

208.0 


frame base strap 

4 

162 

648 

0.104 

67.4 

0.057 

3.8 

13.1 

50.4 


Total Weight (lbs) 



Total weight (lbs) 




878 

878 

projected Area 


48852 

339.25 

sq ft 


1 b/ft A 2 

2.59 



Table 3. Component weight breakdown of a test bulkhead FE model with vertical frames. 


V. Integrated Multibay and Bulkhead 

The Mixed model of integrated multi-bay with bulkhead model, shown in Fig. 8, was studied for two purposes; 
1) Simplify the model to determine an equivalent plate to represent the bulkhead and cargo -bay for 2P load; and 2) 
Capture the critical are area behavior for asymmetric bending -torsion loads. The upper crown panels and side 1 
bulkheads are modeled with foam-core frames and rod stringers. The bulkhead side 2 is modeled with 1.5 inch thick 
equivalent flat plate. The bulkhead side 2 is hidden to show the interior details. Bulkhead side 2 upper and lower 
bays and cargo bays were also modeled with semi-equivalent plate of 1.5 inch thickness. The equivalent thicknesses 
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were selected to provide same order of deflection on both bulkhead side 1 (stiffened thin shell) and side 2 (uniform 
shell). Figure 9 shows the corresponding displacements and principal strain with 18.4 psi internal cabin pressure. 
Although the equivalent plates cannot be used for buckling analysis, they are necessary for scaled system study, 
comparison with theory, and for building simplified full vehicle model for flight load analysis. This model is 
presently being developed. Integrated PRSEUS fuselage and bulkhead models will also be developed farther for in- 
plane load and buckling analysis. 



Rod-stiffener spacing 6 in. 


modeled as beams 


Frames 
spacing 
X 20 inch 


Rod-stiffener 


modeled as 


beam at 8 in 


spacing 


Frame with 
foam-core, 
and 2 stack 
wrap and 
flange- 
modeled 
as beams 


1.5 inch equiv. 
t skin 

Rod-stiffener spacing 6 in. 
-modeled as beams 


Figure 8. Integrated multi-bay fuselage section and pressure bulkhead configuration. 



Disp. with 1.5 inch equiv. shell thickness 


Displacement (inch) Frame, rod- 
stringers (modeled as beam). Side 
walls, passenger floor, cargo bay 1 .5 
inch uniform shell). (Bay4C) 


Displacement 


(inch) 



0.0004 

■ -0.0004 

■ -0.0012 
1.5 inch equiv. shell thickness 


Principal Strain: Frame, rod-stringers 
(modeled as beam). Side walls, 
passenger floor, cargo bay 1 .5 inch 
uniform shell). (Bay4C) 


Principal Strain 
El (inch/inch) 

0.0087 

■ 0.0079 

■ 0.0071 


Figure 9. Integrated multi-bay fuselage section displacement and principal strain with 18.4 psi cabin pressure. 


VI. Summary of FEM Weight Analysis 

The first and second-generation 800 passenger HWB concepts were developed at Boeing Phantom Works 
through NASA contract and NASA in-house research 1 ' 5 . In Ref. 5, four structural concepts were studied, namely, 
Vaulted Ribbed Shell (VLRS), Flat Ribbed Shell (FRS); Vaulted shell with Light Honeycomb Core (VLHC) and 
Heavy Honeycomb Core (VHHC); and Flat sandwich shell with Light and Heavy Honeycomb Core (FLHC & 
FHHC). A multi-bubble section concept was analyzed in Ref. 6 and a relative comparison of component non- 
optimal weight breakdown per unit surface area is shown in the first bar chart in Figure 10. Conventional two floor 
cylindrical and A3 80 type elliptic section and stress balanced multi -bubble fuselage sections with frame and 
longiron stiffened aluminum skin of 3 mm (0.118 inch) shell thickness were also analyzed for comparison. 

The weight analysis summary of the Multibay (without pressure bulkhead) is shown in the 3 rd bar charts in Figure 
10. The weight analysis includes three design scenarios: 1)2 stack skin and rod-spacing of 8 in.; 2) 2 stack skin and 
rod-spacing of 6 in.; and 3) 1 stack skin and rod-spacing of 8 inches. Passenger floor has 1 stack skin and rod- 
spacing of 8 in. The frames are at 20 inch spacing for all 3 cases. The Multibay configuration weight with 2 stack 
skin and 2 stack frame wrap appears be a safe-design for this internal pressure load case. The weight can be reduced 
by using a single stack skin configuration, but the impact damage can be a serious concern. The estimated weight of 
the multi-bay with 1 stack skin and two stack frame wrap is approximately 1800 Ih. (about 1.92 Ib/sq. ft of surface 
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area) for one stack skin and two stack frame wrap, 8 inch rod-spacing). The maximum deflections are of the order 
1.37 inches at outer skin surface of the cargo area (Figure 6). The maximum stresses are of the order 69,000 psi on 
the skin at the corner junction. The maximum strains are of the order 0.007 in./in., but local strains are higher at 
corner junctions. In the 4 th bar chart, estimated specific weights of three bulkhead configurations are presented from 
preliminary analysis of three configurations: 1)2 stack skin, 2 stack rods at 8 in spacing (Fig.7b); 2) 2 stack skin, 1 
stack rods at 8 in spacing (Fig. 7a); and 3) 1 stack skin, 2 stack rods at 8 in spacing. In all case the vertical frames are 
at 24 inch spacing. In the first case, the bulkhead weight is 878 lb (2.6 lb/sq. ft. surface area) without considering 
edge reinforcements and attachment fittings. The maximum defection is about 1 inch. The stresses and strains are 
generally within the allowable limits. Although the weights can be reduced in the 2 nd and 3 rd configuration, the 
deflection and stresses can be significantly higher. 



VHHC+BH+SW FHHC+BH+SW VRBS+BH+SW FRBS+BH+SW 


fi; 

tziTI In 11 

□ 

Qm; 





Cylinder Elliptic 2-Bubble 3-Bubble 4-Bubble 5-Bubble 
7.5m dia. 7x8m 7.5x11. 6m 7.5x15.5m 7.5x19.3m 7.5x17.2m 


O" 

(/) 



PRSEUS Multi-Bay7D 
Wt (IbJJarea (sq ft) 


o- 

(/) 


3.00 

2.50 


1.50 

1.00 


1 


1 . 2 stack skin, rods @8, 5 frames @20 

2. 2 stack skin, rods @6, 5 frames @20 

3. 1 stack skin, rods @8, 5 frames @20 





PRSEUS Bulkhead2T 
Wt (lb) /area (sq ft) 























1 2 3 


1 . 2 stack skin, 2 stack rod@8, frames@24 

2. 2 stack skin, 1 stack rod@8, frames@24 

3. 1 stack skin, 2 stack rod@8, frames@24 


Figure 10. Summary of relative non-optimal finite element model concept weight per unit surface area. 


VII. Conclusions 

This finite element analysis and trade studies of a hybrid-wing-body fuselage section indicate that 
advanced, stitched, composite stiffened plates with rod-stringers and foam-core frames are better than the vaulted 
double-skin-shell or multi-bubble stiffened-shell concepts. A combination of skin-thicknesses, rod-stringer spacing 
and frame spacing were considered for the weight reduction trade study. The finite-element-model based specific 
weights of an ideal three-bay section and bulkhead model are estimated to be 2.0 lb./sq. ft. and 2.6 lb./sq. ft., 
respectively, depending on the configuration. These specific weight values are comparable to that of the cylindrical 
fuselage, but are lower than multi-bubble concepts, for the 18.4 psi pressure load condition. Although, 0.052 inch (1 
stack) outer skin is preferable from a minimum weight perspective, additional load conditions and buckling studies 
indicate that 0.104 inch (2 stacks) skin will provide better safety margins and impact damage resistance. Since the 
sandwich frames take most of the bending load, two stacks of frame -wraps and frame caps are necessary. The 
combined maximum local bending and axial stresses in the rod-stringers also appear to be high at end-joint areas. 
With continuing development and testing, this stitched composite technology could make the current hybrid-wing- 
body concept both structurally feasible and aerodynamically efficient. But, the adverse effect of pillowing on the 
outer-surface aerodynamic flow should be considered. Systems analysis of the stitched composite bi-axially 
stiffened shell concept demonstrated that the hybrid-wing-body class of vehicles can be structurally as efficient as 
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the conventional cylindrical skin- stringer-frame construction. However, alternate concepts should also be considered 
to improve the present design farther. 
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